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Abstract 

Given a set S of n points in general position, we consider all fc-th order 
Voronoi diagrams on S, for k = 1, . . . ,n, simultaneously. We deduce symmetry 
relations for the number of faces, number of vertices and number of circles of 
certain orders. These symmetry relations are Independent of the position of 
the sites in S. As a consequence we show that the reduced Euler characteristic 
of the poset of faces equals zero whenever n odd. 

1 Notation. 



MS) 

B{a,b) 
CH{S) 
SU 
VkiS) 

öa,b,c 
öa,b.c 

[n] 

U{S) 

fi 
Ii 

6; 
Ci 



the arrangement defined by ali bisectors in S. 

the bisector of a and b. 

the convex huU of S. 

the boundary of the set U. 

fc-th order Voronoi diagram on S 

Circle defined by points a, b and c 

number of points from S inside the circle. 

{l,...,n} 

The Voronoi poset on S. 
number of regions in Vi{S). 
number of i-diniensional faces. 
number of vertices in Vi{S). 
number of edges in Vi{s). 
number of circles of order i. 



2 Introduction. 

The dynamics of Voronoi diagrams in the plane is well understood. When n — 1 
points are fixed and one point is moving continuously somewhere inside the convex 
huU, combinatorial changes of the Voronoi diagram correspond to changes in the 



configuration of empty circles. See for example |AGMR|. Changes in the configu- 
ration of non-empty circles correspond with combinatorial changes of higher order 
Voronoi diagrams. Here the fc-th order Voronoi diagram associates to each subset 
of size k of generating sites that area in the plane that consist of points closest to 
these k sites. 

We consider all fc-th order Voronoi diagrams simultaneously for k between 1 and 
n. We do so by introducing the Voronoi poset of a set 5 of n different sites in the 

*This research is financially supportcd by the NWO-stichting SWON, under project number 
613-02-204 
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plane. The poset consists of ali sets of labels that correspond with a subset of sites 
that defines some non-empty Voronoi region in some fc-th order Voronoi diagram. 

Higher order Voronoi diagrams have been investigated by numerous people. Many 



results are published in an article by D.T. Lee, see |Le|. A survey is given in 
Edelsbrunners book on algorithms in combinatorial geometry, see | Ed ] . He ends his 
paragraph on the complexity of higher order Voronoi diagrams with the following 
remark. 

Interestingly enough, the number of regions of Vk{S) is thus exactly 
fcn — fc^ + 1 if n is odd, if fc = {n + l)/2 and if no three sites in S are 
collinear and no four sites are cocircular. 

We generalize this result by giving a symmctry relation for the number of regions 
f k in the fc-th order Voronoi diagram. Assuming general position as in the quote, 
the following cquation holds for any k. 

fk + fn-k+i = 2fc(n-fc + l) + l-n 



See Lemma 5.1. We prove a similar result for the number of vertices, see Lemma 5.2. 
As mentioned above, circle configurations have a close relation with higher order 
Voronoi diagrams. Suppose we are given a set S" of n points in general positions. 
Every three points of S define an unique circle, see Figure ^. Let Ci denote the 
number of circles defined by S that contain exactly i points of S. Then Theorem 



5.2 proves 

Ci + Cn-i~3 = 2(i + l)(n - z - 2) 




Figure 1: An invariant for circle configurations. 

Note that these two symmetry relations are Independent of the particular position 
of the sites in S: although the number of, for example, regions in some fc-th order 
Voronoi diagram will vary, depending on the configuration, the sum of the number 
of regions in the fc-th order diagram and the number of regions in the (n — fc + l)-th 
diagram will stay the same. 

Closely related is the question between which bounds the number of regions of the 
fc-th order diagram can vary. Tight lower bounds for fc — 1 < n/2 are known, see 



Ed|. The symmetry relation gives some bound on the other side. 
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Another application application of the symmetry relations is the analysis of the 
reduced Euler characteristic of the Voronoi poset. It turns out that 

mis)) = o 



whenever n is odd, see Theorem 




As there is a tight connection between higher order Voronoi diagrams in the plane 
and level sets of hyperplane arrangements corresponding with the projection of the 
Voronoi diagram on the unit paraboloid in M^, some of these resuhs may help in 
the investigation of complexity issues with respect to those level sets. 

I most certainly want to acknowledge Wilberd van der Kallen and Dirk Siersma for 
asking me many critical questions guiding me in good directions. 

3 Higher order Voronoi diagrams. 

Definition of fc-th order Voronoi diagram. 

Suppose we are given a set S = {si, . . . , s„} of n points in the plane. Essential in 
this work is the following assumption, the general position assumption. 

general position. 

1. No more than three points froni S lie on a common circle. 

2. No more than two points from S lie on a common line. 

Let O < fc < n. For every point p in the plane we can ask for the k nearest points 
from S. That is, we look for a subset A C S, such that 

1^1 = fc, Va; G A, Vye5-A: d{p,x) < d{p,y) 

For two points in R^, we define a halfplane 

h{x,y) {peR^ \d{x,p)<d{y,p)} 
We define the Voronoi region of A C S" of order \A\ as the intersection of halfplanes 

ViA) := fl hix,y) 

x£A, y£S-A 

whenever this intersection is not empty. As an intersection of halfplanes, V{A) 
will be a convex polygon. Assuming general position implies that a Voronoi region 
cannot degenerate to a line segment or a point. This is proved in Appendix A. 

We define the k-th order Voronoi diagram as the subdivision of R^, induced by the 
set of Voronoi regions of order k. For later purposes, we identify the fc-th order 
Voronoi diagram with the set of non empty fc-th order Voronoi regions. 

VkiS) := {ViA) lAdS, \A\ = fc, V{A) ^ 0} 
The Voronoi poset. 

In the following fix a labelling of the sites in S and identify a set of sites A C S that 
defines a non-empty Voronoi region V{A) with the set of labels L{A) C [n] of the 
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sites in A. Thus, a subset L of [n] might or might not correspond to some Voronoi 
region V{Al). 

Note that for k = 1 we get back the ordinary Voronoi diagram, which means that 
we have the correspondence 

V^iS) - {{1},{2},...,M} 

We define that ^0(5") = {0}. Vn{S) corresponds to the set {{1, . . . , n}}. We consider 
the set of ali Voronoi regions that appear for a given set S of points and call the 
set of corresponding labels the Voronoi poset 11(5') of S. 

U{S) \J{L{A)\V{A)eV,iS)} 

k 

This definition also makes sense when we drop the general position assumption. 



Circles and higher order Voronoi diagrams. 

We State some elementary properties of higher order Voronoi diagrams. 

Every edge in Vk{S) is part of some bisector B{a,b), a, 6 e 5, and the Voronoi 
vertices are exactly those points that are in the circle centres of three points from S. 
Therefore, under our general position assumption, every Voronoi vertex has valency 
three. The following theorem describes the local situation around a Voronoi vertex. 



Theorem 3.1 [De] Let x he the centre of 0^ ^, for a,b,c ^ S. Let 

H = { z € S \ d{x, z) < d{x, a) } 

and let k = \H\. Then x is a Voronoi vertex ofVk+i{S) and Vfe-|_2('S'). The Voronoi 
edges and regions that contain x are given in Figure |^. 

b a 



V(H U {a,b)) 




Voronoi edge of Vk+i (S) 



V(H U |c I) Voronoi edge of Vk+2 (S) 

I 
I 

bla 



Figure 2: The Voronoi diagram around x. 



In fact, ali vertices, edges and regions can be described as in Theorem 3.1, see [De| 
for more details. Let a, 6,c and H be as in Theorem 3.1. 
of a Circle Qa,b,c ^ W\- Notation: | O^ ,, | 



We will define the order 
\H\. Thus the order of a Voronoi 
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Circle ^ ^ equals the number of points of S" — {a, b, c} it contains. Denote the 
number of circles of order k by c^-and the number of vertices in a k-th order Voronoi 
diagram by Vk ■ As a consequence of Theorem |3.1| we get 

Vk = Ck-l + Ck^2 (3.1) 



Counting vertices, edges and regions. 

The foUowing theorem shows that the total number of vertices, edges and Voronoi 
regions does not depend on the positions of the points in S, assuming general 
position. 

Theorem 3.2 Let Vk, e-k and fk denote the number of vertices, edges and regions 
in Vk{S) for some set S of size n in general position. Then the total number of 
vertices, edges and regions in the Voronoi diagram of ali orders are as folloujs. 

(^)- ELi^fc ^ in(n-l)(n-2) 
(ii)- ELi^fe = \n{n- 1)2 
(iii). ELi fk = \n{n^ + 5) 
Proof. 

(i) . Every circle centre defined by three different sites from S* is a Voronoi vertex in 
some fc-th and (fc + l)-th order Voronoi diagram. As there are (g) different circles, 
the first claim follows. 

(ii) . Consider the arrangement of bisectors ^(S*). Fix one bisector _B(a, 6). Without 
loss of generality we can assume that bisector -B(a, b) is dividcd into n — 1 line 
segments by the Voronoi circle centres 

06x3, abxji, . . . , abxn 

where S = {a, 6, X3, . . . , a;„}. Every line segment is an edge in some fc-th order 
Voronoi diagram. As there are (2) different bisectors, claim (ii) follows. 

(iii) . The Euler formula holds for every order. 

Vk - ek + fk = 1 



Thus 



k=l 



Vk 



k=l 



k=l 



This completes the theorem. ■ 

The number of vertices, edges and regions in Vk{S) is not independent of the con- 
figuration of S as already the ordinary Voronoi diagram shows. But the following 
theorem gives expressions for those numbers, depending on n, k and numbers of 
unbounded regions. Let f^ denote the number of unbounded regions in the fc-th 
order Voronoi diagram. By definition f^ := 0. 



Theorem 3.3 Q Let S be in general position. Then the number of vertices. 



edges and regions in the k-th order Voronoi diagram can be expressed as follows. 
(i). Vk = 2{fk - 1) - 
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(il), ek = 3(A - 1) - /~ 
(m), h = (2fc - l)n - {P - 1) - Eli 



Substituting A: = n in the formula for /j. in Theorem 3.3 gives a formula for the 
total number of unbounded regions. Note that /„ = 1. 

n 

^/-i ^ n(n-l) (3.2) 

1=1 

The unbounded regions in the fc-th order Voronoi diagram can be characterized as 
follows. Let Ipq denote ali points on the linc defined by the points p and q that are 
in between p and q. 

Property 3.4 [OBS] A region V{A) of the k-th order Voronoi diagram Vk{S) is 
unbounded if and only if one of the folloming two conditions holds. 

(i) . There exists a line l that separates A from S — A. 

(ii) . There exist two consecutive points p and q, with p, q G S — A, on 6CH{S — A) 

such that the points in A — pq are in the open half plane defined by Ipq opposite 
to CH{S-A). 

As we assume general position, we only have to consider condition (i) 'm Property 
3.4. It is clear that in this case the following symmetry holds. 

/r = fn-k (3.3) 



4 The Voronoi poset. 

Consider the Voronoi poset 11(5') introduced above. We can order the faces in the 
poset by set inclusion of the sets L{A). This gives us indeed a partially ordered set. 
For more on partially ordered sets consult |zj. The poset is bounded as we have 
the empty set as O, the unique minimal element and the set [n] as 1, the unique 
maximal element. In general, a poset is called graded if it is bounded and if every 
maximal chain has equal length. We show that 11(5') is graded. Below we give an 
example that shows that 11(5) is in general not a lattice. 

Property 4.1 11(5) is graded. 

Proof. We show that r(L( A)) = |i(j4)| \s a. rank functionioYl{{S). A rank function 
maps an element x from a poset to a unique level in such a way that the level 
corresponds with the length of any maximal chain from x to 0. Let L{A) G n(5), 
with |L(^)| = k. Then every point x G V{Ä) has the k points from A as its k 
nearest neighbors. Now order those points with respect to their distance to x. As 
we assumed general position it is always possible to change the choice of x in such 
a way that this order is strict. By removing at each step the furthest point still 
available, we get a chain of length k that descends to 0. ■ 

We analyse the two smallest cases, assuming general position. 
n=3. 

For n = 3 we only have one poset, the fuU poset on [3]. That is, 
n3(5) = {0,1,2,3,12,13,23,123} 
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n=4. 



for n = 4 we have two essentially difTerent posets. We can see this by looking at the 
circles defined by the four points. As n = 4, Voronoi circles will have order one or 
two. Because of Eulers formula, we cannot have four circles of order 1. We cannot 
have less than two circles either as this would give us not enough cells in the first 
order diagram. 




Figurc 3: The two different first order Voronoi diagrams. 



Two empty Voronoi circles. 

Look at the left picture in Figure ^. The only subset of [4] that is missing, is clearly 
23, SO we are left with 

n4(S'i) = {0, 1, 2, 3, 4, 12, 13, 14, 24, 34, 123, 124, 134, 234, 1234} 

This example shows that the Voronoi poset is in general not a lattice. For being 
a lattice it is required that every two elements of the poset have a unique minimal 
upper bound. But the elements 2 and 3 have two minimal upper bounds, 123 and 
234. 

Three empty Voronoi circles. 

This situation corresponds with the picture on the right. The region 123 cannot 
appear in the third order diagram, but ali other subsets of [4] do appear, thus 

n4(52) = {0, 1, 2, 3, 4, 12, 13, 14, 23, 24, 34, 124, 134, 234, 1234} 

5 Symmetry relations. 

Given a set S of sites, we can count for every order k the number of vertices, v^, 
the number of edges, Cfc and the number of non empty Voronoi regions, f k- By the 
f-vector of 11(5) we will mean the vector {/i, /2, . . . , /„}■ The c- and e-vector are 
defined analogously. 

Symmetry in the number of regions. 

It turns out that there exists a symmetry in the /-vectors. 
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Lemma 5.1 Consider the f-vector o/ 11(5'), itihere \S\ — n. Then fk + fn-k+i is 
a constant thafs Independent of the position of the points in S . More precisely, 



f k + fn-k+l 



2k{n - fc + l) + l- n 



(5.1) 



Proof. We apply Theorem 3.3 to f k and fn-k+i- 

k 



fk + fn-k+1 = (2Ä-l)n-fc2 + i_^/-^ 

+ (2(n-A: + l)-l)n-(n-fc + l)2 + l- ^ f°^^ 



n-k+l 



k n — k-\-l 

= 2fcn-2fc2 + 2fc + l-n + n(n-l)-(^/~i+ ^ f°°_^) 

i=l i=l 

We join the two sums by applying the Symmetry Equation and evaluate the 
result by using Equation . 

k n—k+1 n 



i=l 



i=l 



The lemma follows from combining the two equations. B 
Symmetry in the number of vertices. 

A similar equation holds for the number of vertices of a collection of Voronoi dia- 
grams Vk{S), for fc = 1, . . . , n — 1. 

Lemma 5.2 Let S he a set of points in general position, \S\ = n. Let Vk denote 
the number of vertices in the k-th order Voronoi diagram. Then 

Vk + Vn-k — 4fc(n ~ k) — 2n (5-2) 



Proof. Using Theorem 3.3 we can write Vk + Vn-k in terms of numbers of faces. 



Vk + Vn-k = 2(A - 1) - /r + 2(/„_fe - 1) - /„°°_fc 

Regroup and apply symmetry, Equation |3.3| . 

- 2(A. + - 2 - /D 



Apply Theorem 3.3 



n — k 



= 2{n' -2n + 2kn - 2^ ~ (E + E /-i + f^^^ 

i=l i=l 

Combine using symmetry again, 

n 

= 2(r? -2n + 2kn - 2^^ - E ft^i) 

i=l 

Use ELi /z-i = ri{n - 1) 

= -2n + Aku - Ak^ 
= 4:k{n - k) -2n 
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Symmetry in the number of Voronoi circles. 



Recall that the order of a Voronoi circle equals the number of points of S it contains 
in its interior. Thus we define the c-vector of S as the vector 

C{S) = {co, Cl, . . . , C„_3} 



where Ci denotes the number of circles of order i. The following theorem states 
that given n arbitrary points in general position, the number of circles that contain 
exactly i points on thcir inside plus the number of circles that contain exactly i 
points on their outside is a constant. We can prove this by applying the above 
results. 

Theorem 5.3 Consider the c-vector o/ 11(5'), nihere \S\ = n. Then Ci + c^^i^z is 
a constant thafs Independent of the position of the points in S . More precisely, 

Cr + Cn-i-z = 2{i + l){n-2-i) (5.3) 
= 2i{n - ?; - 3) + 2(n - 2) 

Proof. We prove the theorem by induction. 

We use inversion. Inversion changes the point-inside-circle relation in 2-dimensional 



Space in a point-below-plane relation in 3-dimensional space. See [BKOS , page 177 



for more details and further references. The inversion map (f) is defined by 

(j>:{x,y) 1-^ {x,y,x'^ + y'^) 

It lifts points in the plane to the unit paraboloid in three-space. As every circle 
defined by S in the plane contains only three points from S, every hyperplane 
defined by 4>{S) contains only three points from 4>{S) as well. cq, the number of 
empty circles of S in 2D equals the number of faces of the lower huU of 4>{S) in 3D 
and Cn-3, the number of circles that contain ali other points of S equals the number 
of faces of the upper huU of 4>{S). Note that ali images of points in S under cf) are 
part of the convex huU of 4>{S). As the convex huU of a p oint set consisting of n 



points consists of 2n — 4 faces, if every face is a triangle, see | BKOS |, Theorem 11.1, 
the claim follows. 

[induction step] 

Cfe + Cn-k-3 = Ck-1 + Cfe + Cn-k-3 + Cn-k-2 — (Cfc-1 + Cn-k-2) 



Apply Equation 3.1 



Vk+l + t'n-(fe+l) ~ (c/c-1 + Cn^k-2) 



Use Lemma 5.2 and the induction hypothesis, 



= 2(2(Ä: + 1) - l){n - (fc + 1)) - 2{k + 1) 

- {2{k - 1 + l)(n - 2 - (fc - 1))) 
= 2{k + l){n - 2 - k) 
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Let fk := fk + fn-k+i and q := q + Cn-i-3- By the reduced f-vector, denoted /, 
we mean the vector of f k' s for ali different k. That is 

/ ■= {/O) /l; • • ■ ; /[ii^J } 



c is defined similaiiy. As a consequence of Lemma 5.1 and Theorem 5.3, / and c 
are only dependent on n. 



Example. 

As an example we present the reduced /- and c-vectors for n S {3,.. .,12}. 



n 


/ 


c 


3 


(4, 6) 


(2) 


4 


(5, 9) 


(4) 


5 


(6, 12, 14) 


(6, 8) 


6 


(7, 15, 19) 


(8, 12) 


7 


(8, 18, 24, 26) 


(10, 16, 18) 


8 


(9, 21, 29, 33) 


(12, 20, 24) 


9 


(10, 24, 34, 40, 42) 


(14, 24, 30, 32) 


10 


(11, 27, 39, 47, 51) 


(16, 28, 36, 40) 


11 


(12, 30, 44, 54, 60, 62) 


(18, 32, 42, 48, 50) 


12 


(13, 33, 49, 61, 69, 73) 


(20, 36, 48, 56, 60) 



Remark. 

Computer calculations did not suggcst any similar symmetry relation for the number 
of edges. 



Relations between regions and circles. 
Corollary 5.4 

fi = /o + Ci-i = Ct-i+n + 1 



(5.4) 



Proof. This follows directly from Lemma 5.1 and Theorem 5.3 



Property 5.5 Let f°° denote the number of unbounded faces in the i-th order di- 
agram and let q denote the number of circles of order i . 



fr + - C.-2) 



2(n - i) 



(5.5) 



Proof. We prove the property by induction. 

[i=l] c_i is zero by definition. The number of vertices vi in the first order Voronoi 
diagram equals the number of circles of order zero, cq. Thus, 

/r + (co"C_i) = /r+«i 



Apply Theorem 3.3 



+ 2(/i - 1) - 
2(n- 1) 
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[induction step] Assume we have proved that 
We can rewrite this, by using induction again, as 



i+l 

1 = 2m-i(i + l)-^^_i (5.6) 
fe=i 



Now we evaluate — Ci_i. 



Cl - Ci_i = + Ci_i) - 2cj_i (5.7) 
= 2(/,+i - 1) - - 2c,_i 



We fill in this expression for q — Ci_i and apply Theorem |3.3| and Equation 5.6 



/^^i + (c.-c._i) - 2(/,+i-l-c,_i) 

= 2(71 - i - 1) 

This proves the claim. ■ 

Corollary 5.6 The c-vector totally determines the f-vector. The correspondence 
is given by 

f k = n - k + l + Ck^2 



Proof. Applying Equation 5.5 we get 



^/-i = (fc-l)(2n-fc)-Cfc_2 



The claim now foUows from evaluating Theorem 3.3 using this expression 



6 Euler characteristic. 

As an apphcation of the symmetry rclations we will investigate the reduced Euler 
characteristic of the Voronoi poset 11(5'). 

By the reduced Euler characteristic we mean the quantity 

n-l 

X(A) := 5] (-1)7. (6.1) 

i=-l 

where denotes the number of j-dimensional faces of the complex or poset A. For 
a polytope P, for example, the famous Euler-Poincare formula states that x(^) = 0. 



Using the symmetry relation 5.1 we can analyse the Euler characteristic of the 
Voronoi poset. Mind the difference in notation. f ^ equals the number of faces of 
dimension i, where fi stands for the number of faces in the i-th Voronoi diagram. 
Thus 

fi ~ fi+l 
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Theorem 6.1 Let S he a set of points in general position, with \S\ — n > 3. 
Assume n is odd. Then the reduced Euler characteristic o/ 11(5') equals zero. 

m{S)) = o (6.2) 

Proof. Write /i = fi+i+fn-i — fi + fn-i-i- Then we get from the definition 
that 

where 

Tl-3 

tn := 

1=1 

f _i = 1, as f_i counts the empty set. /o is the niimber of zero dimensional faces 
plus the number of n — 1 dimensional faces, so fo = n + l. Applying Equation 
gives 

= -(-1) ' —r- 



Straightforward calculations show that 



tn = (-1) — n 



So it follows that 



, .1+1 71^ + 3 , , n+i + 3 
X = -l + n + l-{-l) — — — + {-!)— — n = O 



Remark. 



Theorem 3.2 doesn't hold whcn n is even. But using similar techniqucs onc can 



prove the following. 

n = 0(4) ^ x{S) odd 
n = 2(4) ^ x{S) even 

Note that as Vk = Ck-i + Ck-2 it follows immediately that 

n-l 



fe=i 

for ali n, where Vk denotes the number of vertices in the fc-th order Voronoi diagram. 
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Appendix A. A Voronoi region is non-degenerate. 



Property A.l Assuming general position implies that a Voronoi region cannot be 
a line segment or a single point. 

Proof. Suppose first that V{A) is a line segement. That means that V{A) is 
locally an intersection of at least two halfplanes h{x, y) and h{u, v) such that the 
bisectors B{x, y) and B{u, v) coincide. This impHes that x, y, u and v are cocircular. 
Thafs a contradiction. 

Now assume that V{A) is a single point p. Then there exist at least three halfplanes 
hl, /i2 and hs such that /ii n /12 n /13 = p. We consider three cases. 

(i). hl = hi{a, b), h2 = h2{u, v) and a ^ u,b ^ v. 

Thus a,u G A and b,v & S — A. As V{A) is a single point, it must hold that 



As we have 



{B(a,b)nB{u,v)) C (h(u,b) Ci h{a,v)) 



p e B{a, b) <^ d{p, a) = d{p, b) 
p £ B{u,v) d{p, u) = d{p, v) 



and 



p G h{u, b) <^ d{p, u) < d{p, b) 
p G h{a, v) d{p, a) < d{p, v) 

it follows that 

diPj a) ^ d{p, v) = d{p, u) < d{p, b) = d{p, a) 
But this implies again that a, b, u and v are cocircular. Again a contradiction. 

(ii) , hl = hi{a, b), /12 = h2{a, v) and b^v. 

This implies that a,b and v are cocircular on a circle with centre B{a,b) fl 
B{a,v). We consider hz{x,y). 

• X ^ a and y ^b ^ u. Sce case 1. 

• X = a and y ^b^ u implies that y and a, 6, v are cocircular. 

• X ^ a and y = b implies that x and a,b,v are cocircular. 

(iii) , hl = hi{a, b), h2 — h2{u, b) and a ^ u. 

This implies that a, b and u are cocircular on a circle with centre B{a, b) n 
B{b,u). We consider hs{x,y). 

• X ^ u and y ^b. See case 1. 

• X = a and y ^ b. See case 1. 

• X ^ u and y = b implies that x and a, b, u are cocircular. 
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